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A third quantization formalisni is applied to a simplified multiverse scenario. A well defined quan- 
tum state of the multiverse is obtained which agrees with standard boundary condition proposals. 
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I. INTRODUCTION 

It is usually claimed that the quantum state of the 
universe is described by a wavefunction [li] which should 
account for all the physical information which can be ex- 
tracted from the universe. Two main approaches have 
been followed in order to obtain such a wavefunction of 
the universe. First, it can be obtained as the solution 
of the Wheeler-De Witt equation 0. This can actu- 
ally be seen as a Schrodinger equation in which no time 
derivative appears in order to satisfy the time invari- 
ance required by the general relativity theory [3|. Sec- 
ondly, a path integral approach can also be taken [l|. 
Then, the wavefunction of the universe is given by a sum 
over all geometries and field configurations which can be 
matched with a given value of the field configuration, 
defined in a spatial section of the whole space-time man- 
ifold. We make a Wick rotation to Euclidean time in 
order to obtain well-defined integrals, and then rotate 
back to Lorentzian time to obtain the final results. 

In both cases, some boundary conditions have to be im- 

Sised [l|, [j| • The Hartle- Hawking no-boundary proposal 
represents a universe which is created from "noth- 
ing" , meaning by that the absence of space, time and 
matter. Vilenkin [5[ considered a universe which is also 
created from nothing, but through a quantum tunneling 
transition instead. Vilenkin argued that the idea of a 
universe being created from nothing is not crazier than 
the creation of particle in other quantum theories. Thus, 
it seems appropriate to consider that the universe may 
be created as a quantum fluctuation of the gravitational 
vacuum. Once a universe nucleates in the space-time 
foam [a, 0| J it nis-y bubble and eventually jump into an 
inflationary period '5'| , which is supposed to be the origin 
of our current universe. Therefore, we must consider a 
multiverse in which changes in the topology of space-time 
are allowed. 

A variety of multiverse hypotheses have been recentl 
considered from different cosmological viewpoints 
The speciflc meaning which is ascribed to a single uni- 
verse depends on the formalism of the relevant the- 
ory. Some well-known examples include the Everett's 
many- world interpretation of the quantum theory |9|, 



the chaotic inflationary multiverse (lO| or the landscape 
in string theory 11 1| . among others (for an exhaustive 
review, see Ref. |12l | and references therein). In this 
paper we only consider the case of a multiverse which 
is described by a set of quantum oscillators, each one 
representing a causally disconnected region of the whole 
space-time. 

Topology changes were first claimed to appear in the 
quantum physics of black hole evaporation [13], and so 
first arose as a result of trying to take quantum physics 
seriously as a description of the whole universe. On pure 
cosmological grounds, the transition from a matter dom- 
inated universe into a space-time filled with phantom en- 
ergy provides us with another example of a bifurcating 
topology originated from the big rip singularity which 
splits space-time into two disconnected regions. 

Creation of universes and therefore topology changes 
are naturally contemplated in the third quantization for- 
malism [14]. This consists of a further quantization of the 
wavefunction of the universe similarly to how quantum 
field theory is constructed from the Schrodinger wave- 
function of matter fields. The computations are difficult 
to perform in the general case. In this paper, however, a 
simplified model will be presented in which such compu- 
tations can in fact be carried out. 

Moreover, the general scheme of the third quantization 
formalism can be applied to both, a multiverse made up 
of parent universes and a space-time foam filled with vir- 
tual baby universes }14| . Parent universes are defined to 
be large space-time regions with a Hubble length of the 
same order as the Hubble length of our universe. Baby 
universes are considered to be virtual ffuctuations of the 
metric in the vacuum of gravity, and their contribution 
to calculations of the gravitational field is extremely im- 
portant at the Planck length. The gravitational vacuum 
would be also populated with virtual Lorentzian and Eu- 
clidean wormholes. The former can be seen as solutions 
of the Einstein's equations with at least two asymptoti- 
cally flat regions, connecting two separate parts either of 
the same universe or of two different universes. The lat- 
ter can be considered as Euclidean sectors of a Friedman 
space-time [ij, [l^, whose quantum states can be seen 
as the exponentially decaying versions of the oscillatory 



universes from which they were Wick rotated. 

Furthermore, the vacuum state of gravity might also 
be observable as far as it could induce a loss of quantum 
coherence in the matter fields. The debate was centered 
(see e.g. Refs. [la] and [l3l) around the assumption of 
taking doubly or even multiply connected wormholes in- 
stead of simply connected ones in the quantum vacuum 
of gravity. Equivalently, it can be placed in terms of 
whether virtual baby universes are created as single uni- 
verses or in pairs. In the theory of the quantum multi- 
verse, in which it has been shown that general topological 
changes may occur, the loss of quantum coherence in the 
matter field sector seems to be unavoidable. 

Our aim in this paper is then to apply the third quanti- 
zation formalism given in Ref. |14| to a simplified model 
of the universe, which nevertheless retains the fundamen- 
tal features of the quantum theory when it is applied to 
the universe as a whole. It will be shown that in such 
a model a well-defined quantum state for the multiverse 
can be obtained, which satisfies the usual boundary con- 
ditions. Squeezed states [18], which are usually inter- 
preted as quantum states without any classical analog 
[ISl [l9l |. are found in the context of the multiverse, and 
in particular they appear in the context of accelerated 
universes. This seems to give support to the idea that the 
acceleration is due to distinctively quantum mechanical 
effects |20l.[21|. Furthermore, the third quantization for- 
malism used in this paper shows some advantages respect 
to other approaches which are usually taken in quan- 
tum cosmology. For instance, it can be demonstrated 
that the quantum state of the multiverse may be given 
in terms of the states of a quantum harmonic oscillator. 
Dealing with frequencies instead of potentials simplifies 
the computations. Moreover, neither a perturbative nor 
a semiclassical approximation needs to be taken in the 
model presented in this paper, and therefore the state of 
the multiverse which is obtained can describe the global 
state of the multiverse as a whole. 

Before going any further, however, some caveats and 
comments should be made about the ideas and results 
considered in this paper. First of all, dealing with 
squeezed states in the context of the quantum state of 
the multiverse would, at first sight, seem somewhat ex- 
traneous. Actually, squeezed states can be readily inter- 
preted both in quantum optics and in a space-time foam 
made up of baby universes, which are both deployed in a 
common space-time. Baby universes can in fact be taken 
as tiny particles representing small perturbations of the 
space-time field, whose quantum state affects the vacuum 
state of the matter field sector, and so this theory may 
turn out to be testable. 

It is more difficult to imagine how we could test a the- 
ory of a multiverse made up of parent universes. It looks 
quite different even though wormhole communications 
channels may exist through which microscopic particles 
could travel from one universe to other [22| . Further- 
more, some correlations among the quantum state of dif- 
ferent universes could also be considered. For instance, if 



the creation of universes would be produced in entangled 
pairs, whose correlations might induce some observable 
effects on the state of each individual universe in the pair 

Squeezed states have already been studied in other 
cosmological contexts such as the inflationary universe 
and gravitational waves [2J, [2^ . It was shown [24i| that 
the gravitational vacuum is populated by gravitons which 
evolve into a highly squeezed state. In this paper it is also 
demonstrated in the context of a third quantization for- 
malism that a space-time foam made up of baby universes 
is in a squeezed state. If such a result could be extrap- 
olated to parent universes, these universes might not be 
independent but quantum mechanically correlated, too. 

The paper can then be outlined as follows: A discus- 
sion with the main features of the third quantization for- 
malism has been also added as Sec. II. In Sec. Ill, it 
is obtained the wavefunction for a quantum multiverse 
made up of Friedmann space-times filled with an homo- 
geneous and isotropic fiuid, and the appropriate Fock 
space is defined. We then analyze the two interesting 
limits of the quantum state of both a large parent uni- 
verse and the quantum gravitational vacuum, which turn 
out to be described by a squeezed state with no classical 
analog. In Sec. IV, the quantum state of the multiverse 
is represented by a density matrix rather than by a wave- 
function, accounting thus for mixed states as well as pure 
ones. Three representations are employed: first, the sec- 
ond quantized wavefunction is used as the configuration 
variable, the usual boundary conditions of Refs. [l|, |j| 
are applied and a probability interpretation can then be 
used. Secondly, a general squeezed number representa- 
tion is taken in order to study the quantum state of the 
multiverse. Parent universes will be represented by num- 
ber states. We show that the quantum state for the vac- 
uum of baby universes is given by a squeezed state, and 
that high order correlations appear among them. For the 
sake of completeness, we have also used the usual P{a) 
representation in terms of coherent states. In Sec. V, we 
compare the results obtained in this paper with previous 
works which have considered or criticized the third quan- 
tization formalism. The conclusions are collected in Sec. 
VI, together with some comments to extend the model to 
a two-dimensional wave equation which would explicitly 
account for other matter fields. 



II. THE THIRD QUANTIZATION FORMALISM 

The basic idea of the third quantization formalism is 
14| : to treat the many- universe system as a quantum 
field theory on superspace. However, the name 'third 
quantization' may be misleading. Actually, the proce- 
dure is formally similar to that of the 'second quantiza- 
tion' used in the quantum field theory of matter fields 
(see. Fig. [1]). The name third quantization comes from 
the fact that the field which is quantized is the wave func- 
tion of the universe, and it depends on already second 



quantized matter fields, i.e. it depends on the particles 
existing within the universe. 

Broadly speaking, the third quantization procedure 
starts from the Wheeler-De Witt equation, it takes one 
of the fields in that equation as a time variable and then 
we can proceed as it is usually done in quantum field the- 
ory. However, it not clear at all that such a 'time variable 
field' can be found in the general 3+1 dimensional theory. 
Such a field can be found in one-dimensional cosmolog- 
ical models as well as in 3-1-1 dimensional models with 
high symmetry. Therefore, let us restrict our attention 
in this section to those cases. 
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rather described by a set of points. On the other hand, 
time dependent matter fields within the universes are 
contemplated. Therefore, let us start with the second 
quantized action for such a system. It reads. 



^^ / 1 
dr \ -TrGu 



.Xt^X" -Afr 



(1) 



where X^^, fi — l,...,D represent D matter fields, G^i/ 
is the supermetric in that superspace and Af is the lapse 
function. For simplicity, a mass term has been chosen 
to be the potential of the fields in the action ([T]). Then, 
we can proceed as usual by defining the conjugate mo- 
menta of the fields, P^ = jjG^^X^ , and the Hamiltonian 
density of the system, i.e. 



H = G^'-'PnPu 



(2) 



In the canonical quantization procedure, P^ — )• —ih ^Z^^ , 
and the classical Hamiltonian constraint, H = 0, turns 
out to be, H(f> = 0. This is the Wheeler-De Witt equation 
for the model being considered and then, cf) = (f>{X^^), rep- 
resents the (second quantized) wave function of a single 
universe. 

To third quantize this theory, let us find an action 
whose variational principle leads to the Wheeler-De Witt 
equation. This can be done with the following third 
quantized action. 



FIG. 1: Analogy 
tion schemes 114 



between the second and the third quantiza- 



A. One dimensional universes 

Let us start then with one-dimensional model (1 tem- 
poral -I- spatial dimensions). Throughout this section, 
we are closely following Ref. [ij]. In that model, a uni- 
verse is represented by a point in space. However, if 
topology changes are allowed, a many-universe system is 



= 1 /"d^XV^(G^''V^0V,0 + m202)^ (3) 

where G = detG^^. Variation of the action ^ with re- 
spect to 4> leads directly to the Wheeler-De Witt equation 
and, therefore, all the information of the second quan- 
tized theory is preserved in the third quantization for- 
malism [141]. 

The usual formalism of quantum field theory can 
then be directly applied to the third quantization many- 
universe system. For instance, the three point Euclidean 
Green function can be defined as [iJI, 



J 



Ge{X^,X^,X^) ^ ^X J d^'X^VG Ge(X^,X^)Ge{X^,X^)Ge{X^,X^). 



(4) 



r 



A is the coupling constant and Ge{X^,Xq) is the two tion (j3|) for this interaction scheme, i.e. 
point Euclidean Green function, i.e. 



GE{Xt,Xf) - / Vcj, 0(Xf)0(Xpe 



MA^-Se 



(5) 



Se^\J d^X^ VG (g^'^V^^SI, 



+ to2(/)2-H-(/)3 



(6) 
S'e[0] is the Euclidean version of the third quantized ac- Thus, the third quantization formalism provides us with 



a system of interacting universes. For instance, the prop- 
agator ^ would physically represent the bifurcation of 
one universe (see, Fig. [2]). 

We can also define a wave function for the whole mul- 
tiverse, i.e. 'i'[(j){X'^), X'~'], where X'^ is the field chosen as 
a time variable and then, i — 1, ..., D ~ 1. It is given by 
the solution of the third quantized Schrodinger equation. 



m)=^f^^\^), 



(7) 



where H is the third quantized Hamiltonian (not to be 
confused with the Wheeler-De Witt operator ^). This 
is constructed from the action (jSj with the momentum 
conjugate to the wave function of the universe defined by 



6C 



(8) 



where, 



d4> 



(in the general covariant form, P^ = 
-G G^^Vj/i/)). The meaning of the wave func- 
tion of the multiverse given by Eq. ([7]) can then be un- 
derstood as follows: let, {|n)}|„=i^2,...: be an orthonormal 
basis or number states, then 



|vl/>=^*(X°)|r 



(9) 



where, $(X°) = \['(X°, 0(X*)), can be seen as the prob- 
ability of n universes at time X^ or the probability of n 
universes with field values X^ |14| . 

This formalism for one-dimensional universes can be 
adapted for 3+1 dimensional universes [IJ]. The formal- 
ism has to be changed to account for the gauge symme- 
tries which are not present in the one-dimensional case. 
However, some of the concepts of the one-dimensional 
case can be applied too to the case of more realistic uni- 
verses. Thus, from the point of view of the third quan- 
tization scheme, our single universe is just one universe 
within the whole many-universe system. Nevertheless, in 
some approximation we have to recover the single uni- 
verse we inhabit. This consists of a large region of the 
space-time of order of our Hubble length called a parent 
universe. Nevertheless, within a single parent universe 
quantum fiuctuations of the metric can also be consid- 
ered. Let us notice that the gravitational action depends 



on the length scale and at Planck length, such fluctua- 
tions become of the same order of the metric [g]. Thus, 
tiny regions of the space-time of order of the Planck scale 
can virtually tunnel out of their parent space-time. These 
are called baby universes. 

The parent universe is influenced by such quantum 
fiuctuations of the gravitational vacuum, for instance 
these can be seen as the source of a non-zero vacuum 
energy or cosmological constant, A. In the third quanti- 
zation scheme, therefore, we may say that a parent uni- 
verse propagates then in a plasma of baby universes. It 
can be depicted as follows: let, (j)p and 0f,, be the wave 
functions of a parent and a baby universe, respectively. 
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FIG. 2: Single and doubled lines represent, respectively, the 
propagators for baby and parent universes: a) the bifurcation 
of a baby universe; b) the nucleation of a baby universe in the 
parent space-time. 



The second quantized action for each universe is given 
by Eq. ([T]) with masses mp and niin respectively. If we 
consider that the difference between the mass scale of a 
parent and a baby universes is large, i.e. rUp » nib, 
then, quantum transitions from baby to parent universes 
are exponentially suppressed. However, Feynman dia- 
grams like those represented in Fig. [5] can still be consid- 
ered. The third quantized action which represents then 
a single parent universe propagating in a plasma of baby 
universes is given by. 
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(10) 



Baby universes are not observable by themselves but 
their effects may be observable, so it seems to be bet- 
ter to consider an hybrid action. To do this we deal with 



<f>p as a wave function in the second quantized formal- 
ism and with (p^ as an operator in the third quantized 
formalism. Then we can interpret baby universes as tiny 



particles deployed in the parent space-time. The result is 
an effective second quantized action for the parent uni- 
verse in which a potential term appears accounting for 
the effects of the plasma of baby universes. In such a 
scheme, it can be assumed an effective interaction action 
of the form, 



Si ^ j dTMY,L,{T,: 



"&: 



(11) 



where the index i labels the different modes of the baby 
universe field (i.e. it labels different species of baby uni- 
verses) , and Ci (r, x) is the insertion operator at the nu- 
cleation event. It defines the space-time points of the 
parent universe in which the baby universes effectively 
nucleate. 



B. Homogeneous and isotropic universes 

Let us now sketch the third quantization formalism in 
a more realistic model of the universe. Let us consider 
then the case of an homogenous and isotropic space time 
with some matter fields defined upon it. For the sake 
of simplicity let us assume just one scalar field, ip (the 
generalization to account for more matter fields can easily 
be done). In that case, the second quantized action can 
be written as [261 



S = 



dt[^gAB{q)q^q''-UU{q)], (12) 



where, M is again the lapse function and, q — {cL{t)^ (ys}, 
being a{t) the scale factor. The minisupermetric is given 
then by. 



GAB{a) = 



-4a 
a^ 



(13) 



with, a^ — AttG — ^, being G and m„ the gravita- 

tional constant and the Planck mass, respectively. The 
potential term in Eq. (|12p is given by. 



U{q) = - 



1 
2^ 



+ ^ka' + a^V{^), (14) 



where, k, — —1, 0, 1, stands for hyperbolic, flat and closed 
spatial sections of the space-time, respectively, A is the 
cosmological constant and V{tf) is the potential of the 
matter field Lp. Thus, for homogeneous and isotropic 
models it is clear that the scale factor can play the role of 
a time variable. This can also be seen from the signature 
of the minisupermetric (IT^ . Therefore, the third quan- 
tization procedure can formally be applied in a similar 
way to the case of one-dimensional universes. Then the 
canonical momenta and Hamiltonian, respectively, are 
given by 



PA 
H 



QABq 



dL 

dqA 

^f[\Q^''pAPB+u{q)]^^fH. 



1 

Jf 



(15) 
(16) 



The Hamiltonian constraint, ^ — 0, turns out to be 

T-Lip = 0, where (j) = 0(a, V') is the wave function of the 
universe. Finally, the third quantized action leading to 
the Wheeler-De Witt equation reads. 



^5 = i / da Qg^^ Va0Vb(/. + Ucl? 



(17) 



As a particular example, let us consider a fiuid with 
equation of state given by, p = wp, being w a constant 
parameter, and p and p the pressure and the energy den- 
sity of the fluid, respectively. These are given by. 



1 
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P 
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2N'' 



:V^ + V{p), 
:P^-V{p). 



(18) 
(19) 



On the other hand, for an homogeneous and isotropic 
space-time the equation of conservation of the cosmic en- 
ergy reads. 



dp= ~i{p + p) 



da 



(20) 



It can be easily solved with the help of the equation of 
state, p = wp. The energy density turns out to be then, 
p{a) = caa~^'^^~^^\ Therefore, combining Eqs. p^ and 
P6I) the Hamiltonian constraint turns out to be. 



pl-^a' + na'-p,a-^^+^=Q, 



(21) 



which is the Friedmann equation, being pa = ah. Under 
canonical quantization, pa — >■ ~^^-§^: ^-nd the Wheeler- 
De Witt equation can be written as. 



(/) + a;2(a)0 = O, 



(22) 



where, = 0(a), is the wave function of the universe. Eq. 
((22t can be seen as the classical equation of motion of a 
harmonic oscillator. It can be derived from the following 
third quantized action. 



's^l 



da 



(02 _ ^2(^)^2^^ (23) 



which formally corresponds to the action of a harmonic 
oscillator, where the scale factor a{t) formally plays the 
role of time. The time- dependent frequency, w(a), is 
given by 



uj{a) 



1 Ac 



Kc^a'^ + pqg 



-3ii; + l 



(24) 



where the velocity of light has been introduced for com- 
pleteness. 

In the general case the computations are actually more 
complicated. However, the Wheeler-De Witt equation 
can be seen like a wave equation and therefore, some of 
the developments used in quantum optics might also be 



applied to the third quantized formahsm of the multi- 
verse. For instance, let us make the change, a = In a. 
Then, the Wheeler-De Witt equation can be written as 






(25) 



This can be seen as the wave equation for a radiation 
field propagating in the void (see, Fig. [T|), at the speed 
c — I, where the field is now the wave function of the 
universe, cj) = (j>{a,(p). The source term '21] in Eq. (1251) . 
J{a,(p), is given by 



Jia,v) 



„4a 



ft2 



Viip) 



^2q 



(26) 



Let us take as the boundary condition for the initial state 
of the multiverse that for a zero value of the scale factor 
there is no space, time and matter (i.e. there is 'noth- 
ing'). Analogously, this is what Strominger calls the void 
Il4|, i.e. |0). Then, the state of the multiverse described 



by Eq. (|25|) would evolve as the scale factor is growing 
up into a coherent state given by. 



i^(«))=n" 



Abl-/3fcfcfc 



\0)k. 



(34) 



This is an example of how far can be formally taken 
the analogy between the electromagnetic field and the 
third quantization formalism of the quantum multiverse, 
at least in the case of cosmological models of high sym- 
metry. 



Now, let us expand the modes of the free wave function 
of the universe {J = 0, in Eq. (^5]) ') as. 



Hoi,^) =^bke~ 



-ZOJfcCK-f ifcip 



(27) 



Then, the third quantization procedure eventually con- 
sists of promoting the wave function of the universe into 
an operator, (j). Therefore, let 5j. and bk be, respectively, 
the creation and annihilation operators of universes in 
the mode fc, i.e. 



(28) 
(29) 



Thus, the system described by Eq. (|25p can formally 
be viewed as the radiation emitted by a classical cur- 
rent distribution given by J. Following the analogy, the 
Hamiltonian that would describe the interaction between 
the field cf) and the current J^ is given by |27j . 




V{a)= dif J{a,(p)(j){a,<f). 



(30) 



In the interaction picture, therefore, the state of the mul- 
tiverse evolves (as the scale factor is growing up) into a 
coherent state defined by. 



d_ 
da 



|M//(a))=e-^/°^("')^">/(0)), 



where [27 



e 1 



|/°V(a')da' ^ 



H' 
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/3fe(a) 



da 



d<^ J(a,(p)e'"'="-''='^. 
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(32) 



(33) 



III. THE WAVE FUNCTION OF THE 
MULTIVERSE 

Let us consider a quantum multiverse made up of ho- 
mogeneous and isotropic universes filled up with a per- 
fect fluid, for which the Friedmann equation is given by 
Eq. (|21l) . Following the third quantization procedure, 
the third quantized action for such a many-universe sys- 
tem is given by Eq. (|23|) . It formally corresponds to the 
classical action of a harmonic oscillator, where it is the 
scale factor which plays the role of time. Then, the quan- 
tum state of the multiverse, |^), can be defined as the 
solution of the third quantized Schrodinger equation. 



R\^i,)^ih^\^), 



where 



H 



L 2 



w2(a) 



(35) 



(36) 



Ptj, being the momentum conjugate to the wave function 
of a single universe, (/>. The frequency uj{a) is given by 
Eq. dMl) 

Different disconnected regions of the whole space-time 
may undergo different expansion (or contraction) rates 
because they can be filled up with different energy-matter 
contents and have different spatial geometries, both fea- 
tures being incorporated in the value of the frequency 
uj{a). Therefore, given an orthonormal basis of solutions 
of Eq. ([55]) . ^'^{(p), the general quantum state of a mul- 
tiverse made up of disconnected regions, each being dom- 
inated by a particular kind of energy- matter content, can 
be expressed in terms of a linear combination of product 
states of the form 



*. 



n\i 



)mi 



* 



N„ 



rl 



(37) 



where Ni is the number of universes of type i, represented 
by a second quantized wave function 4>i, which are filled 
up with an energy-matter content described by the fre- 
quency uji{a). Let us notice that the quantum state given 



by Eq. ((37|) is a well-defined state for closed geometries 
and open geometries. 

The frequency uj{a) of the third quantized harmonic 
oscillator which would represent the quantum state of 
the multiverse, contains not just the information of the 
second quantized theory but also that of the classical 
one. Using the Friedmann equation we can show that the 
frequency determines the expansion rate of a single uni- 
verse. For instance, let us first consider the case of a fiat 
de Sitter universe with an exponential expansion. The 
solutions of Eq. (PT|) are then given by, a{t) ~ aoc*""*, 

with ujQ = \/ ^-J^, where the positive branch corresponds 
to the exponential expansion of a flat de Sitter space- 
time. Taking into account the negative branch, the solu- 
tion for a closed de Sitter space-time can be obtained [S] , 
i.e. a{t) = aocoshujot. Furthermore, in conformal time, 
V = f -ITT, the scale factor turns out to be, airi) = -22_ 

' J a{t) ' ' V // cosrj ' 

which can be seen as the two-dimensional metric of a 
cosmological warp drive [21| . Next we consider the anti- 
De Sitter space-time for which A < 0. Analogously 
to the former case, the solution a{t) ~ agcoswoi can 
be obtained, which in conformal time turns out to be, 
o-iv) = — r^ — , and as a final example let us note that 
the Tolman- Hawking closed baby universe used in 21 1 
corresponds to the values, k = 1 and p = 0, in Eq. (1^ . 
In a similar way we can depict other kinds of baby and 
parent universes in the third quantized theory through 
the value of the frequency. 

Quantum mechanically, in the case of a multiverse 
made up of flat universes filled up with a radiation like 
fluid, for which w = -^j the frequency uj{a) turns out 
to be a constant given by w = ^ , and the state of the 
multiverse is represented then by the state of a harmonic 
oscillator with constant frequency which is, of course, 
very well-known. In general, for universes with fluids 
for which w ^ ^, the quantum state of the multiverse 
can be expressed in terms of harmonic oscillators with 
a time-dependent frequency, which have been studied in 
the past [2^, [23] , and to which a renewed attention has 
been recently paid [30[ . Using the method of invariants 
developed by Lewis and Riesenfeld [2^, the solutions of 
Eq. (|35|) can be expressed in terms of the eigenstates 
of an auxiliary operator, /, which is invariant under the 
evolution (i.e. scale factor invariant) given by the Hamil- 
tonian, H, i.e. ^ = *^ff ~ P'^] ~ *^' "which implies 
that, ih-^I\) = H/|), with |) being a Schrodinger state 
vector. Therefore, the action of an invariant operator 
J on a Schrodinger state vector produces another solu- 
tion of the Schrodinger equation. Furthermore, if the 
invariant operator does not involve differentiation under 
the scale factor, relative phases for the eigenstates of / 
can be found in such a way that such states themselves 
satisfy the Schrodinger equation ([BS]) . giving us an or- 
thonormal basis for the space of solutions. Thus, let us 
assume following [28| the existence of a Hermitian invari- 
ant operator, /, with the homogeneous, quadratic form, 
/ = i [ac/)^ -I- /Jp^ -1-7(0,^0}+]. For such an operator to be 



invariant a set of equations has to be satisfied by the coef- 
ficients a, /3 and 7. Solving those equations, the invariant 
operator can be written as, / — i[(-^)0^-|-(i?p0 — i?(/))^], 
where the function R = R{a) obeys the auxiliary equa- 
tion 

R + u'^{a)R--^^Q. (38) 

R 

Eq. (|38p can also be solved in terms of the solutions of the 
Wheeler-De Witt equation, </> -I- w^(a)0 — 0. Let 0i and 
4>2 be two independent solutions of the latter equation, 
then 



R{a) = ^4>l{a)+4>l{a), 



(39) 



is a solution of Eq. ([55]). Moreover, the invariant opera- 
tor / can be expressed in terms of creation and annihila- 
tion operators, such as 



&(«) - v^li + ^^^p^-^-^) 



h\a) 



R 



i{Rp4, - R(t>) 



(40) 



(41) 



which gives, / — h{b^b+ i). For a flat space-time we can 
obtain analytical solutions of Eq. (1551) and R is given by 
Eq. dSH) with 



</)i(a) 




4'2{a) 



/Tra f UJQa'' 
2q 2? 



\ Q^ 



(42) 
where J„ {x) and Yn (x) are Bessel functions of first and 
second kind, respectively, of order n. For other spa- 
tial geometries and energy-matter contents, the Wheeler- 
De Witt equation ([2^ may not be analytically solvable. 
However, it does not introduce any conceptual drawback 
as far as the existence of the solutions of Eq. p2|) is 
guaranteed, and they may be obtained analytically or 
numerically. Nevertheless, analyticity makes clearer the 
computations and for that reason, we will mainly con- 
sider the case of fiat universes for which such solutions 
are given by Eq. ([^^ . In any case the scale factor depen- 
dent creation and annihilation operators given by Eqs. 
(|40| and (|4T|) give the expected results when they are 
applied on the eigenstates of the operator /, \N,a), i.e. 



b{a)\N,a) 

b^ia)\N,a) 

b\a)b{a)\N,a) 



/N\N -1, 
y/N + 1\N- 
N\N,a), 



l,a). 



(43) 
(44) 
(45) 



and, I\N,a) = h{N+^)\N,a). Therefore, the eigenvalues 
are invariant under scale factor transformations. Further- 
more, the scale factor dependent eigenstates, \N,a), can 
be linked to the eigenstates of a harmonic oscillator with 
constant frequency (let us say, uj = 1). In the Schrodinger 
representation, both wave functions are related through 
the following unitary relation, 

^^(0, a) = —i= f/t ^^(^) I , , (46) 



\N,a}, and ^^i^fi) are the eigen- where, 



where, \E'Ar(0, a) s 
functions of the Hamihonian of a harmonic oscilla- 
tor with constant frequency, i.e. Ko^Nif) = f^iN + 
^)^n{'p), with Ho = |p^ + ^<p^. The unitary opera- 
tor, U = U{(j), a), is given by 



U{<j), a) = e 



^ p 2h R^ 



(47) 



with, R = R{a), defined by Eq. 1^. Then, the most 
general solution |^) of the Schrodinger equation (1551) . can 
be written as 29] 

|vl/,a) =^Cwe^""('^)vl/^(0,a)|iV,a), (48) 

N 

where the relative phases, aN{a), are given by 



aN{a)^~{N+-) 



da' 



(49) 



which has a well defined value because the zeros of the 
Bessel functions Jn{x) and Yn{x) do not coincide, and 
therefore R is not degenerate for any value of the scale 
factor. The number eigenstates \N,a) of the invariant 
operator / form an orthonormal basis for the space of 
solutions of the Schrodinger equation, and the eigenfunc- 
tions \I/iv ((?!>, a) in Eq. (^5)) can thus be interpreted as 
the probability amplitude for N universes with a value 
of the scale factor a, for which their second quantized 
wave functions are given by (/). If we consider the quan- 
tum state of a set of universes with different values of the 
scale factor, the probability amplitudes would be given 
by the product of harmonic wave functions (|48|) , for each 
value of the scale factor. However, it seems to be more 
appropriate to consider probability amplitudes regardless 
of the value of the scale factor because, as it has been 
already noticed, the number of universes within the mul- 
tiverse is a scale factor invariant quantity. In that case, 
we integrate over the scale factor to trace it out. On 
the other hand, since (j) has a well-defined semi-classical 
regime [3l|, the quantum state of the multiverse given 
by Eq. ()48|) can represent both a multiverse of homo- 
geneous and isotropic parent universes and a space-time 
foam made up of baby universes. Moreover, let us no- 
tice that unlike in the second quantized theory, in which 
a quantum state of the universe cannot in general be 
defined 32] , the third quantized state of the multiverse 
given by Eq. p8|) is indeed well-defined because there are 
no space-time singularities in the quantum multiverse. 

However, the orthonormal eigenstates \N,a), are not 
eigenstates of the Hamiltonian given by Eq. (p6| for all 
values of the scale factor (except in the case 11) = ^), 
and therefore they are not stationary solutions for the 
quantum state of the multiverse. In fact, in terms of the 
creation and annihilation operators given by Eqs. pO)) 
and (I41L the Hamiltonian (l36l) can be written as 



H = h 



I3_b^+l3+b^ +l3o 



.t.+ ^ 



(50) 



/31 = /3- 



/3o = 




(51) 
(52) 



The quadratic terms in b^ and b in the Hamiltonian given 
by Eq. (j50|) make the quantunLState of the multiverse to 
evolve into a squeezed state [18|. The squeezing effect is 
greater for the quantum state of a multiverse made up of 
accelerated universes [2^ . It might actually be related to 
the quantum nature of accelerated universes [20, [2l| , in 
the sense that squeezed states are usually considered to 
be sharp quantum states with no classical analog [18l.[l9|. 
However, the squeezing effect vanishes as the dominant 
fluid of the single universes that form up the multiverse 
becomes more and more similar to a radiation fluid when 
w —^ ^. Then, the state of the multiverse evolves into 
a set of conventional coherent states, which are usually 
taken to be the most classically possible states in the 
quantum theory. That can also be considered in the case 
of the nuiltiverse, because for w — -^ the wave function 
of a single universe turns out to be a plain wave whose 
argument is the classical action, i.e. 4>rad{a) — e^'S'^=("), 
with Sc{a) = ujoa. Then, the Hamilton- Jacobi equation 
and the classical equation of motion for the scale factor 
are satisfied [3^, being the distribution of momenta a 
delta function centered at the classical momentum, i.e. 
4>rad{p) = 5{p — Pc), with pc = wq the classical momen- 
tum. In that sense, 4>rad represents a classical universe 
and coherent states can associated with a multiverse con- 
sisting of classical universes (see also Ref. |2(j|V 

Therefore, for a general value of the scale factor and 
w 7^ i the state of the multiverse is represented by a 
squeezed state. However, taking into account the asymp- 
totic expansions of Bessel functions for large arguments 
in the case of parent universes with a large value of the 
scale factor, in Eqs. ([5^ and (|42l) . it can be shown that 
R{a) ~ w~2(a), being a; (a) = woa*"^. Thus, the coeffi- 
cients /3± and /3o in the Hamiltonian (I50p can be approx- 
imated to 



n^p- 



i{q~l) 



4a 
/3o -^ uj{a), 



0, 



(53) 
(54) 



so that the Hamiltonian ([50)) turns out to be the usual 
Hamiltonian for a harmonic oscillator in terms of the 
number operator, A^ = b^b, with a scale factor dependent 
frequency given by uj{a). Then, the quantum correlations 
between number states disappear and they turn out to 
be an appropriate representation for the quantum state 
of the multiverse. In particular, at large values of the 
scale factor (i.e.; | 3> 1), the single universe approxi- 
mation for which (TV) = 1, can be considered. On the 
other hand, taking into account the expansion of Bessel 
functions for small arguments in the case of virtual baby 



universes with a small value of the scale factor (typically 

of Planck length) in Eqs. 1^ and (jH]), then, 0i -^ 

_ i _i ■ 

and (/)2 — > J^Q ^ , and therefore, i? — > z^q ^ and i? — > for 

q> 0. Thus, from Eqs. ([51]) and ([52]), it follows that the 

coefficients /3± and /3o can be approximated as 



Pl^/3- 



"4' 



where i^o is a constant parameter given by. 



vq 



7r{2qy 






2ql 



■It) 



(55) 
(56) 

(57) 



if g 7^ 1. For q — 1, the value of i-o actually coincides 
with hr^ujQ, but in that case /3± = in the Hamiltonian 
given by Eq. (1501) and no squeezing effect is present, as 
has been already shown. However, if we assume that 
baby universes are instead vacuum dominated universes 
or closed baby universes, for which w = —1 and w = — g, 
respectively, then the non- vanishing value /3+ and /?_ in 
Eq. ([50|) when the scale factor decreases implies that high 
order correlations between number states are going to 
play an extremely important role in the space-time foam 
of baby universes at the Planck length. It also means 
therefore that the number state representation is not an 
appropriate representation of the quantum gravitational 
vacuum. Instead, it would be more properly represented 
by a squeezed state, which might be even experimentally 
observable because such a gravitational state would in- 
duce a loss of quantum coherence in the vacuum state of 
the matter field sector [TT'l. In fact, from the Hamilto- 
nian (|50p one can clearly infer that baby universes should 
be continuously created and annihilated in pairs, at least 
under the assumptions made in the present model. That 
appears to enhance the result that euclidean doubly con- 
nected wormholes copiously nucleate, together perhaps 
with simply connected wormholes. This implication is 
relevant for the debate on the loss of quantum coherence 
in the matter field sector due to the gravitational vac- 
uum (see Refs. [l3| and [la]) because it suggests that 
the loss of quantum coherence of matter fields might be 
unavoidable. 

On the other hand, the third quantized vacuum en- 
ergy, Ej, — -|2.^ retains a marginal non zero value even 
when the scale factor decreases to a = 0, i.e - when the 
universes are annihilated. This is a well-known feature of 
the quantum vacuum state with no classical analog, ap- 



plied in the present case to the third quantized vacuum 
state of gravity. 

Finally, let us briefly discuss the question on how 
and when the universes are quantum cosmologically cre- 
ated. In the third quantization formalism, an interaction 
scheme can be used to describe the space-time foam of 
a parent universe. As it was shown in Sec. II, an in- 
sertion operator is then introduced to define the space- 
time points of the parent universe in which the creation 
and annihilation operators of baby universes are located. 
Thus, the interaction Hamiltonian can in general be writ- 
ten as, Hint = J2i fii^T^)9iK^^i)^ where the index i ac- 
counts for the different species of baby universes being 
considered in the space-time foam. 

However, in the case of a multiverse made up of parent 
universes, which cannot be placed in a common space- 
time, the above question becomes misleading. Let us 
notice that space and time are properties with a well- 
defined meaning only within a single parent universe, and 
they are not therefore well-defined outside the universe. 
So the Hilbert space corresponding to a multiverse made 
up of parent universes cannot be described in terms of 
space-time coordinates but, instead, of statistical ones. 
However, correlated pairs of universes can still be con- 
sidered because entanglement is a quantum property di- 
rectly stemming from the superposition principle, and 
therefore, a pair of parent universes can be created at 
points of such a general statistical Hilbert space. 



IV. THE DENSITY MATRIX 

A. The second quantized wave function 
representation 

The most general quantum state of the multiverse is 
not represented by a pure state defined by a third quan- 
tized wave function, but rather by a mixed state de- 
scribed in terms of a density matrix operator, p, for which 
the explicit scale factor dependence is given in the third 
quantization formalism by the Schrodinger like equation, 
-^{p) = jr{[H,p]). Therefore, let us consider the quan- 
tum state of a multiverse made up of flat space-times 
in the representation defined by the value of the second 
quantized wave function </>, i.e. the coordinate represen- 
tation. Then, taking into account the expression of the 
Hamiltonian given by Eq. p6|) . the Schrodinger equa- 
tion for the density matrix in such a representation can 
be written as, 



J 



dp{(j),(t)',a) _ i_ 
da h 



ft2 

2 



a^ 



92 



d(jP d(j)'^ 



r 



w2(a) 



'Z''^) I P{<P,<P' ,0)- 



(58) 



This can be solved by using the following Gaussian 
ansatz, 



p{(j),(j)' ,a) 



-A(a)(^+4>'f-iB(a)(4>^-4,'^)-C(a) 



(59) 



Inserting Eq. (|59p into the Schrodinger equation (l58l) , we 
obtain the following equations for the coefficients, 

A = -AhBA, (60) 

C = 2hB, (61) 



B = -2hB' —. 



(62) 
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The last of these equations can actuahy be transformed 
into the Wheeler-De Witt equation by making the follow- 
ing change: X = e'^'^f^'^" (i.e. B = ^f ). Then, Eq. ^ 
turns out to be x + u'^x = 0, the solutions of which are 
given by Eq. (|42p . The quantum state of the multiverse 
can therefore be expressed in terms of the two values of 
the second quantized wave functions of the universe, cf) 
and (/)' respectively, as 



p{(j),(j)\a) = -— e 

00 



K0+0')' 



i 't'O , 

"2R *n > 



(63) 



where the constants of integration have been chosen 
to satisfy the normalization criteria, Jd(/)|p((?!), 0, a)| = 
1, Va, and, ^o = 4>oifl) = -^o4>i{o) + ^o02(a)i is the gen- 
eral solution of the Wheeler-De Witt equation. Aq and 
B{j are two constants to be determined by the boundary 
conditions. The trace of the density matrix given by Eq. 
(j63| shows no divergences even when 0o — >■ 0, it vanishes 
for an infinite value of the scale factor, avoiding in this 
way cosmic singularities, and it also satisfies the usual 
boundary conditions imposed to the second quantized 
wave function of the universe when the scale factor de- 
generates as a — T' 0. This can be checked by taking into 
account the limits of the Bessel functions in Eq. ((42]) 
for a vanishing value of the scale factor, i.e. 4>i{a) -^ 

and 02(a) ~^ ^o^ ■ Thus, if Eq = 0, the solution of the 
Wheeler-De Witt equation is given by 0i (a), and the den- 
sity matrix (|63|) tends to zero as the scale factor vanishes, 
satisfying in this way a boundary condition which corre- 
sponds to the Hartle-Hawking's no-boundary condition 
[l[. On the other hand, if i?o 7^ 0, then 0o — ^ Bqv~2 ^ and 
the density matrix (|55)) still represents a small oscillation 
even when the scale factor decreases to a = 0, satisfy- 
ing then a boundary condition which corresponds to the 
Vilenkin's proposal of a tunneling boundary condition 
U . The state of the multiverse is therefore well-defined 
by Eq. (|M|) . and the diagonal elements, |p(0, 0, a)|, can 
be then interpreted as a probability distribution which, 
for a given value 4>, is peaked around the solution of the 
Wheeler-De Witt equation, 0o(a)j for which the distribu- 
tion |p(0, 0, a)| is a maximum, i.e. 90op(0, 0)10^=0 = 0. 

Let us recall that the state of the multiverse repre- 
sented by the density matrix (j63p avoids cosmic singular- 
ities, at a = and a — 00^ because it gives a zero value of 
the probability amplitude both when the scale factor de- 
generates and when it grows up to infinity. For instance, 
let us consider the case of a phantom energy dominated 
universe in which a bip rip singularity appears in the fu- 
ture of our proper time, splitting the universe into two 
causally disconnected parts before and after the singu- 
larity. In the context of the third quantized multiverse, 
such a model is rather described by two universes, pro- 
vided that their quantum states are well defined by Eq. 
(|63l). which might be even entangled |3l|, and therefore 
no singularities appear in such a multiverse. That is an 
example of the way in which the third quantized for- 
malism can also account for topological changes in the 



space-time, like the one that would happen, for instance, 
in the transition from a matter dominated universe into 
a phantom dominated one (i.e. from one single universe 
to a pair of, may be entangled, universes). 

Let us also notice that the quantum state represented 
by the density matrix (|63p is a non-factorizable state for 
any value of the scale factor. However, it cannot be in- 
ferred that the quantum state of the multiverse is then in 
a mixed state. It would only be so if the representation 
given by the second quantized wave function of a single 
universe is orthogonal. Nevertheless, the solutions of the 
Wheeler-De Witt equation are not in general orthogonal 
|3l| , and therefore it is not at all clear in such a represen- 
tation whether the multiverse stays in a pure or a mixed 
state. For instance, for large values of the scale factor the 
quantum state represented by the density matrix (p5)) is 
still a non-factorizable state. However, it has been shown 
that in the same limit the number state representation 
is orthogonal and diagonalizes the Hamiltonian. Thence, 
the quantum state of a single parent universe is described 
by a pure state defined by the value of the quantum num- 
ber, N = 1. The number state representation could be 
seen then as an orthogonal or decoherenced coarse grain- 
ing [3JI, given by a superposition of second quantized 
wave functions weighted along accordingly to the third 
quantized wave functions of a harmonic oscillator. 

Finally, let us recall that the density matrix for the 
quantum state of the universe was already obtained in 
Ref. [331 • Our consideration of the density matrix in this 
section can be just viewed then for the aim of complete- 
ness concerning the properties of the quantum state of 
the multiverse in the third quantization formalism. 



B. The squeezed number representation 

Let us consider now the number state representation 
defined by the relations given in Eqs. (|43ll45p . It has 
been shown that in the case of parent universes with a 
large value of the scale factor, or equivalently in the semi- 
classical limit where f — >■ cxi, such a number representa- 
tion diagonalizes the Hamiltonian H, because ;9± — > 
in Eq. ((50)) . In that case, the scale factor dependence 
of the density matrix, p — J^n m -^NM\N,a){M,a\, is 
given by the Schrodinger equation which implies a dif- 
ferential equation for the density matrix elements, i.e. 
iPNM = Po{M — N)Pj\fM, where the dot stands for the 
derivative with respect to the scale factor. The solutions 
are given then by. 



VM(a) = e-'(^-^)/"'*'^'^"(''VwM(0), (64) 



P, 



where /3o(a) is again defined by Eq. ((5^. with /3o(a) ~ 
uj{a) = ^^^^ — for large values of the scale factor. There- 
fore, if the initial density matrix is diagonal in such a 
number representation, i.e. Pnm{0) = 6nmPn, then 
the diagonal elements remain constant with respect to 
the scale factor. For instance, for the pure state given 
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by Pn = Sin, which would represent a large single uni- 
verse, there are no off-diagonal elements in the represen- 
tation of the density matrix as the scale factor continues 
to increase and, therefore, once a large universe can be 
considered, transitions to other number of universes are 
asymptotically suppressed along the subsequent expan- 
sion of such a single universe. 

However, for small values of the scale factor in such 
a representation, non-diagonal elements appear in the 
Schrodinger equation and the number states turn out 
to be squeezed states. Let us note that in terms of 
the creation and annihilation operators for modes of a 
unity value of the frequency, 5g = -^{4> ~ ip<t>) and 
60 = ^^(0 + iPtp), the creation and annihilation opera- 
tors given by Eqs. (^0]) and (jH]) can be written as, 

b{a) ^ ^060 + >^bl, (65) 

b\a) = ^Ibl + v^bo, (66) 



where. 



Mo 



vq 



2 \r 



iR 



R-iR 



(67) 
(68) 



with, l^oP — Wo\^ = li and R = R{a) defined by Eq. 
([55)1 . It means therefore that the number eigenstates of 
the invariant operator / do not diagonalize the Hamilto- 
nian of the multiverse, and the squeezing effect becomes 
then larger for small values of the scale factor. However, 
let us also notice that quantum effects like squeezing or 
entanglement crucially depend on the chosen representa- 
tion [36| . This can be easily shown in the realm of quan- 
tum optics with a simple example: one photon state in a 
given mode of constant frequency, namely wi , which is of 
course a number eigenstate in the above representation, 
turns out to be a squeezed state in the representation of 
the number states given for another different frequency, 
say UJ2, i.e. 61 = /i52 + i^t»2' with |/xp — |i^p = 1, being 61^2 
and b[ 2 the annihilation and creation operators of each 
mode, respectively. In quantum optics the appropriate 
representation is derived from the experimental setting, 
e.g. if the experimental devices of the setting measure the 
number of photons with energy Huid, then the eigenstates 
of the number operator of that mode, N^ = bdbd, would 
become the appropriate representation. The choice of the 
representation is not that clear nevertheless in the case 
of the quantum multiverse. After all, identifying the rel- 
evant physical observables in the multiverse is not at all 
an easy task. In the case of the quantum multiverse, 
it can be assumed that the Lewis number states, which 
are defined as the eigenstates of the invariant operator /, 
represent the occupation number for the quantum state 
of the multiverse, i.e. the number of universes in the 
nrultiverse, which is scale factor invariant. 

However, other representations can be more useful in 
order to obtain the solutions of the Schrodinger equa- 



tion for the elements of the density matrix, i.e. ^ (p) = 
jr{[H,p]). For instance, let us consider the diagonal rep- 
resentation of the Hamiltonian (|36| . It can be defined 
by the scale factor dependent creation and annihilation 
operators which correspond to the modes of the Hamil- 
tonian with proper frequency, w(a), i.e. 



buj{a) 
blia) 



uj{a) 

2h 



w(a) 

i 
UJ(a) 



(69) 

(70) 



where, uj(a) 



!£0„q-l 



In such a representation, of 



course, the Hamiltonian turns out to be diagonal and 
given by, H = tiLu{a){bljb^ + i)- The annihilation and 
creation operators given by Eqs. ([69)1 and ([70)) can, more- 
over, be related to the annihilation and creation opera- 
tors of modes of unity frequency, bo and 6j, through the 
following squeezing transformation. 



buj = Mtj^o + t^Ljbl, 
bl = f^lbl + iy^bo, 



where, 



Mw 



i^,., = - 



(71) 
(72) 



(73) 
(74) 



and hence, fi^ 



1. Then, the new number states. 



\N^, a), can be related to the number states of modes of 
unity frequency, \No), by 

\N^,a)^S^ie^)\No), (75) 

where the squeezing operator, Scj{Si^), is defined by 

5'„(£^) = e^^"(^o-''o'). (76) 

Ei^ is the squeezing parameter which is related with Eqs. 
([751) and (TMl) by, ^^^ = coshe^^, and, v^ = sinhe^. As 
a matter of completeness, let us note that the number 
eigenstates of the Hamiltonian can also be related to the 
Lewis number states, \N,a), defined in Sec. I by Eqs. 
([ ^ - ([l5)) . as they are both linked to the number states 
of modes of unity frequency, \No). Thus, ii Sii{eii) is the 
squeezing operator that relates the Lewis number states, 
\N,a), with those of unity frequency, |A*o), defined by 



SBiSR) 






(77) 



^2i0 _ 



where, eji — re^*"^, with, r ~ cosh/ig, and, o — -. — r, 

being /io and vq given by Eqs. ([67)1 and ([68)1 . then 
\N,a) = Sji\Nq). Therefore, the eigenstates of the 
Hamiltonian, \N^), are related to the number states. 
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I TV, a), through a more comphcated squeezing transfor- 
mation given by, 



\N^,a)=SU£cu)Sl^{eR)\N,a} 



(78) 



The squeezing transformations are unitary so that 
S^s) = S~^{e) = S{—e), and therefore, they preserve 
the orthonormahty relations among the number states of 
a given mode. Thus, the Hamihonian eigenstates, \Ni^), 
form also an orthonormal basis for the space of solutions. 
In order for these to satisfy the Schrodinger equation, 
they have to be shifted by some phases which can be ob- 
tained following the same iterative, recurrent procedure 
as the one used in Ref. |28| . 

Then, the orthogonal representation of Hamiltonian 
eigenstates that satisfy the Schrodinger equation is given 
by, 



|^^,a)=e*(^+iM-)|iV^,a), 



where, 



n{a) 



Wo 



LO[a )aa — -—a 
nq 



uj{a) 



q 



(79) 



(80) 



Let us now compute the scale factor dependence of the 
density matrix elements by solving the Schrodinger equa- 
tion in such a representation of Hamiltonian eigenstates, 
\Ni^, a), in which the density matrix operator can be writ- 
ten as, p = X^AT A/-PiVMl-^'^'")(-^^'^'^l- The solutions to 
the Schrodinger equation are given then by. 



P^Mia)^e'^('''-^)r 



da 



"^">]^m(0), (81) 



where, a = a; -I- iy, is a complex valued variable, and 
il{a) is given by Eq. (|80|) . The coherent states (|83|) are 
the eigenstates of the annihilation operator for a mode 
of frequency w, with scale factor eigenvalues given by. 



b^{a)\a,a) = e*"(°)a|a, a). 



(84) 



In the customary P{a) representation, the action of the 
creation operator upon the coherent state (|83| is given 

by, 



6i|a,a)=e-''"('') 



_d_ 
da 



a, a), 



(85) 



and the density matrix can be written then as, 

p{a)= f (fa\\a,a){a,a\\e-\"^^P^{a,a), (86) 



where, ||a, a) = e^i~\a,a), are the so called Bargmann 
coherent states [33|. Thus, using Eq. (|55)) in the 
Schrodinger equation, daP — j^[H,p], a differential equa- 
tion is obtained for the function P{a, a), i.e. 

d f d d \ 

—P^{a,a) = -iuj{a) la- ^* JTT j Puj{a,a). (87) 

This equation can be easily solved by making the change, 
a — X -\- iy. The solutions can be written as. 



P^{a,a) =CiPi(a,a)+C2P2(a,a). 



(88) 



Ci and C2 are two constants to be determined by the 
boundary conditions and Pi (a, a) and P2(a, a) are two 
independent solutions of the master equation (|87p , given 
by 



with the diagonal elements taking on constant values. 
By employing another number representation it is possi- 
ble to re-express the density matrix in terms of squeezing 
transformation operators. For instance, in terms of the 
number states of modes of unity frequency, |TVo), the el- 
ements of the density matrix operator are given by. 



P^Aa) ^ (/oIpIJo) 



(82) 



E 

N.M 



P^M («) (-^0 1 S^ (e„ ) I iVo) (Mo I Si (£„ ) I Jo) 



where the involved matrix elements can no longer be triv- 
ially calculated. 



C. Coherent state representation 

Let us consider next the P representation of the quan- 
tum state of the multiverse in terms of coherent states. 
These can be defined in terms of the number eigenstates 
of the Hamiltonian as. 



00 

E- 

7V=0 



..N 



/N\ 



3*(^+5)"(")|iV^,a), (83) 



Pi (a, a) = -(e-*f^('^)a 



e,ma)a* 



P2(a,a) 



1 L-m^^) 



a + e 



in{a) 



(89) 
(90) 



where J7(a) is once again given by Eq. (|80l). and the inte- 
gration constants have been chosen to ensure the real 
valued character of the function P^. Eqs. ((89|) and 
([^O]) closely recall the form of the operator which rep- 
resents the potential vector, AEMix^t), in the usual sec- 
ond quantization formalism of the electromagnetic field 
|38| . Thus, Eq. ([55]) may be interpreted as an oscillating 
field of a particular scale factor dependent mode, ri(a), 
although the quantized field is now the wave function of 
a single universe. 

Other solutions to Eq. ([57)) can be found as well. 
For instance, by expanding in series the function P^, i.e. 
P(j = J^ni^riOi" + Bna*"'), the solution can be written 
as, 



P^{a,a) = Re 



1 



l±e-'^('')a 



P^(0,0), 



where n{a) is given by Eq. (|80|) . and the integration 
constants have been again chosen to ensure a real value 
for the function P,,. 
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In terms of other number representations, it may be- 
come difficult to carry out tlie explicit computation of 
the function P. For instance, let us consider a coherent 
state representation in terms of a general number repre- 
sentation,, i.e. 



p{a) 



where. 



d^a\\ag,a){ag,a\\e-\''\"Pgia,a), (91) 



OO ff 

,a)^5:^e^(^+^)"«('^)|7V„a). (92) 



N=0 



^Nl 



\Ng, a) are the number states of modes of a generic fre- 
quency, ujg = ujg{a), defined in terms of the annihilation 
and creation operators. 



(93) 
(94) 



where for convenience we have retained the (plausi- 
ble) complex character of the frequency of the mode. 
The phases in Eq. ([M]) . which are proportional to 
fig (a), account for the evolution of the number states 
like Schrodinger vectors. Therefore, in such an unspec- 
ified representation the Schrodinger states evolve into 
squeezed states, and the Hamiltonian is transformed into. 




H 



h(j-bl+J+bl'+j0{2blbg 



1 



(95) 



where j± and 70 are functions which in general depend 
on the scale factor, except in the case of the diagonal 
representation of the Hamiltonian. They are given by. 



7- =7+ 



70 



IWoPwo 



4(Re^g)2 



4(Re^g)2 



-f 1 



(96) 
(97) 



and therefore, ^g{a) = J ^^da. It can be checked that 
when the frequency of the representation, Wg, equals 
the proper frequency of the Hamiltonian, w, then such 
a representation diagonalizes the Hamiltonian. In gen- 
eral, the Schrodinger equation for the density matrix, 
daP = Ji[H,p\, can be written as. 



d^a||ag)(ag||e '"' daPgia^a) 



(98) 



- / d2c,[i/,||ag)(ag||]e-l"l'9aPg(a,a).(99) 



Taking into account that the relations ([84]) and (|85l) are 
satisfied now with 5g, w and ^g{a) instead of b^, b\^ and 



^{a), the differential equation for the function Pg{a,a) 
can be written as. 



d_ 
da 



Pg{a,a) 



d l_d d_ 

dxj "' 2 dxi dxj ^■' 



Pg{a,a), 



(100) 

where, (a:i,X2) = (a, a*), the drift vector, Aj, is given 
by 

Aj = {Aa,Aa*) = (C+a* +i7oa,C-a- J7oa*)- (101) 
The diffusion matrix, Dij , is given by 



A,= 



D,, 







i^a-a- 



e+ 
e- 



(102) 



with, ^1 = ^_|_ = 2z7+e~^*^s, being, j'H = 7+. In general, 
Eq. (|100p is not easy to be analytically solved. However, 
before going any further it should be first clarified the 
meaning of the different representations for a particular 
cosmological model in terms of what can be observed, 
which is not at this stage clear at all. 



V. COMPARATIVE ANALYSIS WITH 
RESPECT TO OTHER PREVIOUS WORKS 

The third quantization formalism has already been 
used in previous works on the multiverse and it has been 
criticised (see Ref. [39|, and references therein). First 
of all, Vilenkin argued that an infinite number of fields 
and interaction types had to be introduced in order to 
describe topology changes in a four-dimensional realm, 
i.e. - for all the non-simply connected space-times that 
can generally be considered. 

Let us notice, however, that such a non-simply con- 
nected manifold can actually be converted into a simply- 
connected one by cutting it with a finite number of hy- 
persurfaces J13j|, and this argument can be applied in the 
four-dimensional case. Hawking pointed out [13] that the 
quantum state of the universe would then be given by a 
set of quantum oscillators, an idea similar to the multi- 
verse scenario. In the third quantization formalism, on 
the other hand, the wave function which represents each 
of the disconnected regions might be given by that of a 
harmonic oscillator. In that case, the third quantization 
formalism adds a manageable procedure to study a mul- 
tiverse made of disconnected regions, at least in the case 
of high symmetry. 

The second main objection of Ref. [39| was on the 
interpretation of the creation of universes in the third 
quantization formalism. Vilenkin argued that the usual 
interpretation in terms of \in) and \out) states, corre- 
sponding to \N — 0) and \N — 1), respectively, with 
N the number of universes, does not match the mean- 
ing of the creation of a universe from nothing proposed 
by himsef some years before [40|. The \in) state would 
correspond to the "nothing state" when the scale factor 
decreases, whereas the \out) state would describe a single 
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universe with a large value of the scale factor, i.e. when 
a — >■ oo. However, as Vilenkin himself objected [S^, the 
scale factor cannot be used as a time variable because it 
is not a monotonic variable. Moreover, taking a particu- 
lar value of the scale factor is nothing but considering a 
particular time reparametrization, a — a{t), and it is not 
expected that the number of universes in the whole mul- 
tiverse would depend on a particular choice of a reference 
system within a single universe. 

On the other hand, we have shown that the number 
state which properly defines the number of universes is 
represented by the Lewis number state defined by Eq. 
(|45p and thus, the number of universes in the whole mul- 
tiverse is scale factor invariant, i.e. N ^ N{a). However, 
the concept of universes created from nothing can still be 
meaningfully kept because the vacuum state of the mul- 
tiverse, i.e. the void^ still represents the absence of space, 
time and matter, and is what Vilenkin means when he 
uses the term " nothing" . 

Finally, it was considered in Ref. [41| that in the con- 
text of eternal inflation a mosaic structure made up of 
causally disconnected thermalized regions should appear 
with different values of the effective coupling constants, 
with a probability distribution for their values. However, 
it is not at all clear that the general state of the multiverse 
has to be a thermal state because, unlike the inflationary 
model, a common space-time is not guaranteed at all in 
the case of a general multiverse. 

In our simplifled model, on the other hand, a vanish- 
ing cosmological constant seems probable, in agreement 
with the result of Ref. [43]. However, there may be a 
solution to this problem, because we have shown that 
in terms of the creation and annihilation operators for 
single universes, b^ and b, a non-zero value of the vac- 
uum energy should be ascribed to each universe, even at 
a = 0. Therefore, although the most probable value for 
the cosmological constant would be zero in the realm of 
the whole multiverse, if it is taken to be formed up by 
pairs of universes, a non-zero value for the vacuum energy 
could be present in each single universe due to entangle- 
ment with its partner, which could be responsible both 
for the currently observed accelerating expansion and for 
the value of the vacuum energy in the inflationary era. 



VI. CONCLUSIONS AND FURTHER 
COMMENTS 



be given in terms of a wavefunction or a density matrix, 
accounting for pure and mixed states. Thus, the third 
quantization scheme can be completely well-defined in 
these models of high symmetry. 

The state of the multiverse turns out to be given by a 
squeezed state with no classical analog. On the one hand, 
the quantum state for a gravitational vacuum which is 
filled with a large number of popping baby universes and 
wormholes connecting them, is given by a squeezed state 
in which quantum correlations among the number states 
arc dominant. This result agrees with that of Ref. [2J|, in 
which the vacuum state of gravity is populated with relic 
gravitons which unavoidably evolve also into a highly 
squeezed state having no classical analog. On the other 
hand, in the case of large parent universes the quantum 
correlations among the number states disappear asymp- 
totically, the number state representation turns out to be 
an appropriate one and the single universe approxima- 
tion can be then considered. Thus, once a single parent 
universe can be assumed, quantum transitions to larger 
number states are asymptotically suppressed. Notwith- 
standing, if universes are created in pairs, then the corre- 
lations in the quantum state of a pair of universes would 
imply that each single universe is not independent. Then, 
one of the main non-independence consequences might 
be the presence of an entanglement energy, which could 
account for both the large value of the cosmological con- 
stant in the inflationary period and for the small value 
of the vacuum energy now, because the entanglement be- 
tween the pair of universes decreases as each single uni- 
verse expands. 

Finally, let us notice that the model presented in this 
paper can be naturally extended to include local matter 
fields as well as cosmological fields. In the first case, it 
is expected that their effect would represent a minor cor- 
rection to the value of the frequency of the harmonic os- 
cillator which quantum mechanically described the mul- 
tiverse. However, in the latter case, derivatives with 
respect to the matter fields appear in the Wheeler-De 
Witt equation, implying that the state of the multiverse 
is given in terms of the modes of a wave equation in which 
the mass and the frequency might be scale factor depen- 
dent. Nevertheless, it seems that it will be easier to work 
with the solutions of such a wave equation rather than 
with the general solutions of the second quantized theory. 



In this paper, we have applied a third quantization 
formalism to the case of a multiverse made up of homo- 
geneous and isotropic space-times filled with a perfect 
fluid. A well-deflned quantum state of such a multiverse 
has been obtained in terms of the states of a harmonic 
oscillator, and it satisfles the most familiar boundary con- 
ditions which are used in quantum cosmology. It can also 
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